Josephson effect in S/F/S junctions: spin bandwidth asymmetry vs. Stoner exchange 



Gaetano Annunziata/ Henrik Enoksen,'^ Jacob Linder,'^ Mario Cuoco/ Canio Noce/ Asle Sudb0^ 

^ SFIN-CNR, 1-84084 Fisciano (Salerno), Italy 
Dipartimento di Fisica "E. R. Caianiello" , 
Universita di Salerno, 1-84084 Fisciano (Salerno), Italy and 
Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

(Dated: April 27, 2011) 

We analyze the dc Josephson effect in a ballistic superconductor/ferromagnet/superconductor junction by 
means of the Bogoliubov-de Gennes equations in the quasiclassical Andreev approximation. We consider the 
possibility of ferromagnetism originating from a mass renormalization of carriers of opposite spin, i.e. a spin 
bandwidth asymmetry. We provide a general formula for Andreev levels which is valid for arbitrary interface 
transparency, exchange interaction, and bandwidth asymmetry and analyze the current-phase relation, free en- 
ergy, and critical current, in the short junction regime. We compare the phase diagrams and the critical current 
magnitudes of two identical junctions differing only in the mechanism by which the mid-layer becomes mag- 
netic. We show that a larger number of — tt transitions caused by a change in junction width or polarization 
magnitude is expected when ferromagnetism is driven by spin bandwidth asymmetry compared to Stoner mag- 
netism. Moreover, we show that these features can be present also for ferromagnets of the Stoner type having 
only a partial bandwidth asymmetry. 

PACS numbers: 74.50.+r,74.45.+c,72.25.-b,74.25.Dw 



I. INTRODUCTION 

The antagonistic nature of ferromagnetism and singlet su- 
perconductivity notwithstanding, coexistence of these phe- 
nomena is not entirely ruled out and their interplay may in 
fact yield quite rich physics. Coexistence of ferromagnetism 
and superconductivity, predicted to occur in a phase where 
both the gauge and the continuous translational symmetries 
are spontaneously broken simultaneously,^!^! js suggested to 
emerge in many theoretical models.^"* However, a decisive 
experimental verification of such coexistence is still lacking. 
On the other hand, the interplay between ferromagnetism and 
superconductivity in heterostructures has been the foundation 
for developments of numerous novel experimental techniques 
over the last decade or so.^ ^ Besides the interest from a fun- 
damental physics viewpoint, such systems hold potential in 
terms of future spintronics^ applications. The versatility of 
heterostructures involving superconductors and ferromagnets 
stems from the basic fact that Andreev reflections,'^ ^ processes 
taking place at the interfaces converting quasiparticl es int o 
Cooper pairs, are strongly affected by spin polarization.ESEI] 

In this paper, we focus on ballistic superconduc- 
tor/ferromagnet/superconductor (S/F/S) Josephson junctions 
with conventional superconducting leads. These are systems 
capable of sustaining a supercuiTent caiTied by Cooper pairs in 
the superconducting leads and by quasiparticles in the ferro- 
magnetic mid-layer. The unique interplay between feiTomag- 
netic and superconducting orders provides quasiparticles with 
extra phase shifts absent in junctions with a non-magnetic 
mid-layer. This can give rise to the appearance of a so-called 
"tt— phase". Under such circumstances, the energy mini- 
mum of the junction is reached at a phase difference of tt 
across the junction, unlike the standard "0— phase" in junc- 
tions with a non-magnetic mid-layer. The existence of the 
TT— phase has been experimentally confirmed.Q3Sl As far as 
potential applications are concerned, tt -junctions are consid- 



ered to be candidates for realizing solid state qubits.E^ 

Previous theoretical studies on S/F/S junctionPl22l have 
described the F layer using the Stoner model of metallic fer- 
romagnetism, with oppositely polarized carriers occupying 
rigidly shifted bands. However, the interplay of Coulomb re- 
pulsion and the Pauli principle driving a metal into a ferro- 
magnetic state,!^ may induce a spin-dependent renormaliza- 
tion of the masses of charge carriers with opposite spins, i.e. 
a spin bandwidth asymmetry. This mechanism appears in mi- 
croscopic approaches where off-diagonal terms of Coulomb 
repulsion, generally neglected in studies based on the Hub- 
bard model, are taken into account. A mean-field treatment of 
these contributions in addition to the exchange and nearest- 
neighbor pair hopping terms shows that quasiparticle ener- 
gies for the two spin species are not simply splitted but get 
different bandwidths, i.e. effective masses. The net effect of 
these interactions is to render the hopping integral in the ki- 
netic term spin dependent through bond charge Coulomb re- 
pulsion terms which are different for spin up and down carri- 
ers.l^^For low enough temperature and depending on Hamilto- 
nian parameters, ferromagnetic order can be established only 
through this spin bandwidth asymmetric renormalization.l^Sl 
In this picture the ferromagnetism should be understood as 
kinetically driven, in the sense that it arises from a gain in 
kinetic energy rather than potential energy, unlike the usual 
Stoner scheme. We notice that this mass-split metallic fer- 
romagnetism has been experimentally found to be the ori- 
gin of the optical properties of the colossal ma^netoresistance 
manganites,^^ of some rare-earth hexaboride^ as well as of 
some magnetic semiconductors.'^ Spin bandwidth asymme- 
try will substantially affect the coexistence of ferromagnetism 
and superconductivity,!^ as w ell as the proximity effecP and 
transport in F/S bilayers.'^^'^It is also responsible for an ex- 
tension of the regime in which FFLO (Fulde-FeiTell-Larkin- 
Ovchinnikov) phase can be stabilized in heavy-fermion sys- 
tems.^*-- In this work, we analyze the consequences of this 
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unconventional magnetism in ballistic S/F/S Josephson junc- 
tions by solving the Bogoliubov-de Gennes equations^** for 
arbitrary spin polarization and interface transparency. 

We show that bandwidth asymmetry in the F layer of a 
S/F/S junction modifies many physical properties such as An- 
dreev levels dispersion, Josephson current and free energy. 
We provide a general formula for Andreev levels which holds 
for arbitrary interface transparency, exchange interaction, and 
bandwidth asymmetry and we calculate the corresponding 
Josephson current in the short junction regime for some im- 
portant limiting cases. Comparing the results for a spin band- 
width asymmetry ferromagnet with those obtained for a con- 
ventional Stoner ferromagnet, we demonstrate analytically 
that the former manifests features that are only quantitatively 
different from the latter for low degrees of polarization. We 
show numerically that these differences become qualitative in 
the intermediate/high polarization regime. Then, evaluating 
the phase diagrams we infer that the mere mechanism by itself 
is able to shift the ground state superconducting phase differ- 
ence from to TT or viceversa, with a spin bandwidth asym- 
metry ferromagnet tending to increase the number of possible 
— vr transitions in a fixed range of parameters. 

This paper is organized as follows. In Section II, we intro- 
duce the model and notation, and provide a general formula 
for Andreev levels in S/F/S which is valid for arbitrary inter- 
face transparency, exchange interaction, and bandwidth asym- 
metry. In Section III we present an analytical form for An- 
dreev levels and Josephson current in different limiting cases 
and numerically analyze the same quantities, critical current, 
and free energy for different polarization, temperature, F layer 
width, and interface transparency values. This analysis is per- 
formed considering both Stoner and spin bandwidth asymme- 
try magnetism in the F layer The relationship between the 
spin polarization and the magnetic microscopic parameters is 
clarified by means of a general formula holding for a generic 
ferromagnet in which both mechanisms are present. In this 
"mixed" case, we also show that the features distinguishing 
pure spin bandwidth asymmetry from pure Stoner ferromag- 
nets may exist even though the spin bandwidth asymmetry 
contributes only partially to spin polarization. We present 
phase diagrams and analyze — tt transitions which are driven 
by variations in the polarization and width of the F layer. The 
conclusions and a discussion on the applicability of our results 
to real materials is presented in Section IV, while analytical 
details for some quantities used in the paper are reported in 
the Appendices. 



II. MODEL AND FORMALISM 

We consider a ballistic Josephson junction composed of 
a ferromagnetic layer of width L sandwiched between two 
conventional s— wave singlet superconducting electrodes (see 
Fig.[T]i. Since the S electrodes are isotropic, we will consider 
an effective one dimensional model. Moreover, our results are 



independent on the direction of the magnetic moment direc- 
tion in F. The propagation of quasiparticles is described by 
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FIG. 1: A sketch of the current-biased S/F/S junction analyzed in 
the paper, along with notations that are used. The superconductors 
are treated as reservoirs, whereas the magnetic mid-layer has a finite 
width L. 



the Bogoliubov-de Gennes equationJ^ 



where a — —a, p^(i) = +1(— 1), and (ucVg) = 4'cr is the 
energy eigenstate in the electron-hole space associated with 
the eigenvalue e. The single-particle Hamiltonian is 



m = HL + H^ + Hn + Hi, (2) 



where 

Hl 
HI 
Hr 
Hi 



[-h^\l'^ /2m- Ep] Q[-x) 

[-n^\i^/2m„ - p^u - Ep] e{x)e{L 

[-h'^V^/2m - Ep] e(x - L) 

H [S{x) + Six - L)] . (3) 



x) 



Here, different effective masses for cr-polarized particles 
in the F layer, mimicking bandwidth asymmetry, has been in- 
cluded and U is the exchange interaction, Ep is the Fermi 
energy, Q{x) is the Heaviside step function, m is the effec- 
tive mass of the quasiparticles in the superconductors, and 
the parameter H quantifies scattering strength at S/F and F/S 
interfaces. We assume rigid superconducting order parame- 
ters with equal gap amplitude on both sides of the junction, 
i.e. A = Ao [e"^^e{-x) + e*'^«e(a; - L)] , where (^^(j^) is 
the phase of the left(right) superconductor. When consider- 
ing finite temperature properties, we will assume the usual 



BCS dependence and let Aq Aq tanh (^1.74^Tc/T - 1 
where Tc is the superconducting critical temperature. 

Employing the quasiclassical Andreev approximation,!^ the 
solution of Eqs. ([TJ may be written as 
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^ — ikx 



X < 



X > L 



where 



2e 



2e 



(4) 



and A; = \/2mEp /h, q„ = yj2ma{Ep + PaU) jh are Fermi 
wavevectors in the S and F electrodes respectively, and Lp = 
'PR ~ 'PL- The coefficients a^, b„, c„, and d„ are the probabil- 
ity amphtudes for Andreev reflection as a hole-like quasipar- 
ticle (HLQ), normal reflection as an electron-like quasipar- 
ticle (ELQ), transmission to the right electrode as an ELQ, 
and transmission to the right electrode as an HLQ respec- 
tively, while the coefficients a^, Pa, la, and 5a are associated 
with right- and left-going ELQ and HLQ propagating through 
the ferromagnetic layer. All these probability amplitudes can 
be calculated by solving the linear system resulting from the 
boundary conditions 
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yjUL-)--i''a{L+), (5) 
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which are valid for a general ferromagnetic electrode display- 
ing both exchange splitting and bandwidth asymmetry. When 
the latter mechanism is present, electron- and hole-like parts 
of the wave function derivatives have to be joined separately 
with different masses. This is analogous to the situation where 
the insulating barriers are polarized with magnetic moment 
parallel to the one in the F layer and thus act with differ- 
ent strength on particles with opposite spin.ESlWhen only ex- 
change splitting is present, i.e. nia — rng = rn, the usual 
form of the boundary conditions is recovered.G3 



Inserting the wave function into the boundary conditions, 
we obtain a homogenous system of linear equations for the 
scattering coefficients. By imposing that the system has non- 
trivial solutions, we can obtain a relation between energy 
and phase difference such that the coherent subgap processes 
can effectively take place, i.e. Andreev levels with dispersion 
£i{(p), i = {1, 4}.™We are able to provide a general form 
of the levels which holds for arbitrary transparency of the in- 
sulating barriers, polarization in the ferromagnetic layer, and 
relative weight of the exchange and bandwidth asymmetry 
contributions. We find £2 (1^9) — £a{'P') — ^^cri'P) ^^d 



Eaip) = Ao 



A-^ + B'^ -A{C + D cos2(95/2)) + Pa^B^ l^A^ + B^ - {C + D cos^{ip/2)) 

\ 2 (A2 + B2) 



(6) 



where A,B,C,D depend on all junction parameters. Their 
explicit form will be given in the Appendix |A] Here, we 
merely note that for transparent interfaces, and in the absence 
of exchange field and bandwidth asymmetry, Eq. (|6]l reduces 
to the well known result e — Aq cos{(p/2)F^ 



We will focus on the short junction regime, i.e. L <ti ^ 
where ^ is the superconducting coherence length, such that 
the Josephson current can be estimated by only considering 



subgap contributions by Andreev level: 



as 



li'P) 



dip 



h ^ \ 2 J dp 



(7) 



where f{x) is the Fermi-Dirac distribution function and /? = 
l/keT. 

Finally, the critical current is defined as 1^ ~ max\I{p)\ 
and the phase of the junction (0 or vr) is determined by the 
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minimum of the dependent part of the free energy 
1 



In 



n 1 



2 cosh 



-/3£,(V) 



(8) 



III. RESULTS 

In this section we analyze Andreev levels, current-phase 
relation, free energy, and critical current for a short S/F/S 
junction. Our model enables us to take into account both 
bandwidth asymmetry and exchange interaction in the F layer. 
Defining the polarization as M = (n^ — n-i)/ (n^ + ti^) and 
integrating densities of states as calculated from the Hp term 
in Eq. ([3]), we find that in one dimension and at T = 0: 



M = p+/p^ with Pa 



1 + U/Ep 



1 - U/Ef 



(9) 



For a fixed value of the exchange splitting U, the effect of 
mass mismatch is to enhance the net polarization for to^ > 
mi, and to hinder it the other way around. Eq. (|9| describes a 
general F where both exchange and bandwidth asymmetry are 
present. We notice that Eq.(|9]) reduces to the case where only 
exchange interaction is present when / = 1 and U 
or the case where only bandwith asymmetry is present when 
™t /''^l 1 U — 0. We will refer to these two cases 
as a Stoner ferromagnet (STF) and a spin bandwidth asym- 
metry ferromagnet (SBAF), respectively. When both mecha- 
nisms are present there are several combinations of exchange 
energy and mass ratio that give the same polarization M. We 
choose to describe this mixed case through a parameter W 
e [0, 1] which quantifies the relative weight of the two mech- 
anisms without changing the total polarization in such a way 
that a pure STF and a pure SBAF will correspond to = 
and W — 1, respectively. The detailed procedure chosen 
to calculate the value assumed by exchange and mass mis- 
match for a given polarization M and degree of mixture W 
is illustrated in Appendix |B] In Table |l] the magnetic param- 
eter values for a STF, a SBAF, and a mixture of the two with 
W = 0.5, are reported for three values of the polarization 
M = {0.25,0.50,0.75}. Different Josephson effect features 
are expected in S/STF/S and S/SBAF/S junctions, even for an 
equally polarized F layer This is so, since specifying a degree 
of polarization is equivalent to fix the ratio of Fermi wavevec- 
tors q^/qi- However, in a SBAF and in a STF the wavevectors 
are different, i.e. q^^^ ^ q§^^^, and accordingly the center 
of mass momentum acquired by Cooper pairs will be differ- 
ent. 



In the following, we will also consider finite temperature 
properties in the range < T < Tc- Even if Eq. (|9| was 



M 




U/Ef 






STF 


8/17 


1 


0.25 


SBAF 





25/9 




STF >U SRAF 


1 0^ 

u. i yj 


i .o / o 




STF 


4/5 


1 
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STF & SBAF 


0.287 


4.988 




STF 


24/25 


1 


0.75 


SBAF 





49 




STF & SBAF 


0.324 


24.99 



TABLE I: Values of the normalized exchange interaction and mass 
asymmetry used for three different polarization values. Here STF 
& SBAF represents a mixture of the two magnetic mechanism with 
W = 0.5. 



derived for T 0, it can be used for T < Tc, since the polar- 
ization varies only slowly on the scale of the superconducting 
critical temperature as long as Tm ^ Tc, with Ti\j being the 
Curie temperature. This condition is easily fulfilled in a typ- 
ical experimental situation where the Josephson junction is 
realized by a ferromagnetic transition metals compound sand- 
wiched between conventional superconductors, e.g. Niobium. 
In our analysis, we fix Aq = 1 meV and Ep = 5 eV, consis- 
tent with the Andreev approximation. The width of the ferro- 
magnetic layer L is not fixed. However the maximum value 
used in our analysis is such that L/(_ sa 0.01, such that we 
are in the short junction regime, corresponding to a width in 
the range 1-10 nm. We introduce a dimensionless parameter 
quantifying the scattering strength at insulating interfaces as 
Z = 2mH/ (tfk). 

Before exploring the general case, we look closer at the low 
polarization limit for the Andreev levels to get an idea of what 
to expect. Expressing the mass mismatch in a SBAF and the 
exchange energy in a STF as a function of polarization yields 




1 + Af2 ■ 



(10) 
(11) 



At first glance, Eqs. ( [T0] l and ^H} indicate that exchange in 
a STF only depends on odd powers of the polarization M, 
while mass mismatch in a SBAF depends on all powers. Fur- 
thermore, the expansion of Eq. ( [TOj l is strictly positive, while 
Eq. ([TTJ has alternating signs. Due to the intricate dependence 
on mass mismatch and exchange energy of the general An- 
dreev levels Eq. (|6]l (see the Appendix for details), both SBAF 
mass mismatch and STF exchange will depend on all powers 
of AT However, the different sign behavior may indicate that 
SBAF and STF behave differently for large degrees of polar- 
ization. From now on to lighten our notation we will refer to 
S/SBAF/S and S/STF/S as simply SBAF and STF The first 
evidence of different behavior is encountered in the M ^ 1, 
Z < 1 limit. Expanding Eq. ^ around U/Ep < 1 for STF 
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and yjir^jral 1 for SBAF, we obtain 




(13) 

for SBAF when Z ^ 0. Then, for low polarization and includ- 
ing the first term from the Z <^ 1-expansion, both expressions 



reduce to the same form 

= Ao (cos C^^j + p„M sin [kL + 9Z sin^ kL)^ , 

(14) 

where = +1(-1) for STF (SBAF). Thus in the low polar- 
ization regime we expect SBAF and STF to show the same 
behavior for Z = 0, while the ^-factor indicates that they can 
be different for nearly transparent interfaces. This is a conse- 
quence of a different interplay between the two ferromagnetic 
mechanisms and the insulating barrier strength described by 
Eq. (|5]). By differentiating Eq. ([14]) with respect to if and in- 
serting it into Eq. Q, we obtain the current-phase relation 



^(V') _ 2 *i sin((/j/2) sinh[/3Ao co&{f/2)] - S2kLM cos{(p/2) sinh[^AofcLM sin((^/2)] 
lo ~ cosh(/?£^/2)cosh(/3e;/2) 



(15) 



where si = sgn (cos {(p/2)), S2 = sgn (sin {ip/2)), and Iq = 
eAo / h. We have again put Z = for simplicity. 

From Eq. ( fTSj ), we see that the current-phase relation con- 
tains a cosine term in addition to the standard sine term. One 
might think that the cosine term will create a non-vanishing 
current at ip = 0, tt, but this is not the case. The cosine term 
will cancel out at 1^9 = 0, tt due to the Fermi-Dirac distri- 
bution function, as shown in Eq. ( [T5| ). However, in a non- 
equilibrium situation where levels can be differently popu- 
lated, these terms may give rise to exotic features such as frac- 
tional AC Josephson effect"*^ or a lyS-junction-i^This possi- 
bility is matter of current investigation. 

In the tunneling limit it is possible to obtain an analytical 
form for the current for an arbitrary degree of polarization. 
For T — and including only the first term in the Z ^ 1 
expansion, we find 

I{tp) _ 1 - A/2 2sin(v5) 
^ " 1 + M2 



sin kL 



(l+M)^ 



sin kL 



for STF and 

lo 



2 sm{Lp) 



Z4 sin ( kL 



l + M 
1-M 



sin ( kL 



l-Af 
l + M 



(16) 



(17) 



for SBAF. At first glance, one may think that the main differ- 
ence between the two cases is the term (1 — M'^)/ (1 + AP) 
which is present only for STF and should damp the current 
for high degrees of polarization. We will later see that such an 
effect exists for low Z. However, for high Z this is masked by 
the fact that the current is small. In fact, the main difference 
between Eqs. ([16]) and ( [TT] ) is that the denominator becomes 
wildly oscillating for large M for SBAF, but not for STF. We 
will see later that this effect is clearly manifested also beyond 
the tunneling limit. 

We now move on to describe the S/F/S junction properties 
upon leaving aside the limits of low polarization and high/low 



M=0.25 Z=0 



M=0.25 Z=1 




FIG. 2: (Color online). Positive branches of spin-split dispersion of 
the Andreev levels for different values of polarization and interface 
barrier strength and for Lk — 10, Ef = 5 eV, Ao — 1 meV. Solid 
lines show levels in SBAF and dashed lines in STF. 



interface transparency. Fig. [2] shows the ip dependence of 
and ei in SBAF (solid lines) and STF (dashed lines) for trans- 
parent (left panels) and insulating (right panels) interfaces, for 
three values of polarization M. All levels are flat at ip = 0, tt, 
which corresponds to the absence of Josephson current for 
these phase difference values. This is a well known prop- 
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SBAF M=0.50 




FIG. 3: (Color online). Current-phase relation and free energy for 
STF (left panels) and SBAF (right panels) for intermediate and high 
polarization values. In each panel three different interface barrier 
strenghts are considered: Z — (solid lines), Z — 1 (long-dashed 
Mnes) Z — 2 (short-dashed lines). Here Lk = 10, Ef = 5 eV, 
Ao = 1 meV, T/T,. = 0.01 are fixed. 



erty of S/N/S and S/STF/S junctions, which persists in the 
S/SBAF/S case. However, the dispersion of the levels are 
different for STF and SBAF. This difference is only quantita- 
tive for low polarizations, e.g. M ~ 0.25, as expected from 
Eq. ( [T4] i. However, the levels differ qualitatively for higher 
polarization in terms of their slope and curvature. The order 
relation between and can also be different for the two 
magnetic mechanisms in the F layer, e.g. for M = 0.75 and 
Z = 1 £^ < for STF while > for SBAF 

These differences affect the Josephson current and free en- 
ergy, as shown in Fig. [3] where, I/Iq and F/ Aq are plotted for 
different polarization values, different magnetic mechanisms, 
and Z = (solid lines), Z = 1 (long-dashed lines), Z — 2 
(short-dashed Unes). Here Lk = 10, Ep = 5 eV, Aq = 1 
meV, and T/Tc — 0.01. For low polarization values (not re- 



ported in the figure) the current and the energy differ only 
quantitatively for STF and SBAF, while for intermediate/high 
polarization values a qualitative difference between the effect 
of the two mechanisms emerges. For instance, at M = 0.75 
the Josephson cuiTent exhibits a maximum at < tt for STF 
and 93 > TT for SBAF (see panels c) and d)). This means (as 
shown in panels g) and h)) that the energy minimum of the 
ground-state for the two types of junctions can be realized at 
different superconducting phase differences even though the 
polarization in the F layer is the same. Moreover, we see that 
STF is in a 0— phase while SBAF is in a tt— phase, so that the 
magnetic mechanism by itself is able to drive — tt transitions. 
In the following, we will refer to the situations where SBAF 
and STF are both in 0— or tt— phase as "in phase", while when 
one is in the 0— phase and the other is in the tt— phase as "out 
of phase" junctions. By inspection of panel a), b), e), /) 
of Fig. [3] it is also possible to infer that the interface trans- 
parency can play a role in this bandwidth asymmetry driven 
phase shift: the STF and SBAF junctions are out of phase for 
Z = and become in phase for Z = 1,2, since SBAF under- 
goes a — TT transition for some Z in [0, 1]. We underline that 
STF can also undergo such transitions driven by a variation in 
the strength of the barrier, depending on L and T values. 

In order to ascertain if the realization of out of phase junc- 
tions is a rarity or a common situation, we have compared the 
free energy in SBAF and STF for different width, tempera- 
ture, degree of spin polarization, and interface transparency 
values. Our analysis shows that the condition of out of phase 
junctions is mainly determined by spin polarization: in the 
low polarization regime, i.e. M < 0.25, SBAF and STF junc- 
tions very often share the same phase, while in the intermedi- 
ate/high polarization regime, i.e. AI > 0.50, they are likely 
to be out of phase. Other parameters only weakly affect these 
conditions except interface transparency which favors out of 
phase situation whenever the tunneling limit is considered. 

We now consider the case of a mixed F with both exchange 
energy and bandwidth asymmetry. We show that the previ- 
ously discussed features distinguishing SBAF from STF as 
seen in Fig. [3] are also present when considering a mixed F. 
Indeed the Josephson cuiTent and the free energy for mixed F 
look very similar to the case of a pure SBAF, when a suitable 
contribution of bandwidth asymmetry to total polarization is 
reached. Fig.|4]shows the case of a mixture of SBAF and STF 
with W — 0.5 and all other parameters fixed as in Fig. [3] In 
the mixed case both exchange and masses ratio have non triv- 
ial values. The values assumed by these parameters for the 
particular choice W — 0.5 are reported in Table [I] As previ- 
ously discussed for M — 0.75, STF is in a 0— phase while 
SBAF is in a tt— phase. It is clear from panels b) and d) of 
Fig.|4]that the particular mixed case considered with W = 0.5 
is in a tt— phase just like the pure SBAF case (corresponding 
toW = 1). Depending on the junctions parameters, a mixed F 
can behave as a pure SBAF even at smaller bandwidth asym- 
metry contributions. For the particular case considered here, 
the mixed F behaves as a pure SBAF for W > 0.3. By com- 
paring SBAF and STF at M = 0.5, we have pointed out that 
STF and is in a tt— phase while SBAF is in a 0— phase for 
Z — 0, while for smaller interface transparency SBAF switch 
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FIG. 4: (Color online). Josephson current and free energy for a fer- 
romagnet which is a degree W = 0.5 mixture of STF and SBAF. All 
parameters are fixed at the same values of Fig.|3] 



FIG. 5: (Color online). Critical Josephson current as a function of L, 
M, and T for SBAF (solid lines) and STF (dashed lines). 



to the same phase of STF. Panels a) and c) of Fig. |4] show 
that even in this case a partial bandwidth asymmetry in F is 
enough to observe features of a pure SBAF. For the polar- 
ization considered (M = 0.5), a mixed F behaves as a pure 
SBAF for W > 0.4. These results are a clear signature of the 
fact that the bandwidth asymmetry can induce unusual fea- 
tures also when it is not the only mechanism producing the 
spin polarization. We stress that this is an important consid- 
eration since, as discussed later, in real ferromagnetic mid- 
layers both exchange and bandwidth asymmetry mechanism 
may be present. 

Fig. |5] shows the Josephson critical current as a function of 
L, M, and T for SBAF (solid lines) and STF (dashed Unes) for 
Z = 0,1, 2. When the L dependence is considered with a fi- 
nite polarization, e.g. M = 0.25, the current for both the STF 
and SBAF junction appears oscillatory for transparent bank- 
ers and displays complicated patterns when Z ~ 1. These 
patterns repeat periodically when Ic is observed on a larger 
scale and collapse to a modulated Dirac comb for Z ^ 10. 
The maximum critical currents for the two junctions are very 
close, but they are attained for different L values and appear 
to be weakly affected by the Z value. The M dependence of 
the critical current shows a damped oscillatory behavior for 
SBAF and STF for Z = and more complicated patterns for 
Z — 1,2. The main difference between the two junctions is 
that for high polarization values, AI > 0.6, the critical current 
decreases monotonically for STF, whereas for SBAF it oscil- 
lates rapidly with increasing frequency. This is a signature 



of the fact that SBAF still undergoes — tt transitions, while 
STF has settled in either a 0— or a tt— phase. This behavior 
is consistent with Eqs. ( [T6] l and ( [T7] i even if they are strictly 
valid only in the tunneling limit. The plots of the temperature 
dependence of the ciitical current in Fig. |5] show no partic- 
ular features. The critical currents in both STF and SBAF 
are decaying functions of temperature, but the cuiTent in the 
latter junction is larger. This is not a general feature but de- 
pends on the M and L values chosen. In general one can have 
both I^^^^ > /f^^ and I^^^^ < if^^ as clearly shown by the 
plots of critical current as a function of these parameters in 
Fig. |5] Moreover, a temperature driven — tt transition can 
be seen both for STF and SBAF by carefully tuning M and 
i,''^-e.g. M = 0.675, Lk = 11.4 for SBAF and M = 0.935, 
Lk = 11.4 for STF. Even a slight change in one of the two 
parameters suffices to destroy a temperature driven transition, 
and the current simply decays as in Fig. |5] 

The analysis of the free energy and critical current has 
shown that for a given set of junction parameters such as de- 
gree of polarization, width, and insulating barrier strength, the 
Josephson effect in SBAF and STF can be different. These 
two kinds of ferromagnets can support different critical cur- 
rent magnitudes, can induce different ground state phases 
across the junction, and can undergo a different number of 
— TT transitions driven by variations in junction width and 
polarization. In order to study these points in more detail, 
it is convenient to fix a range for these parameters and then 
look at the behavior of STF and SBAF in this entire range. 
We employ this strategy because the observables are rapidly 
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FIG. 6: (Color online). Density plot of AJ = (jf ^ '^^ - /f^^)//o 
in the parameter range < M < 0.75 (g) 6 < Lfc < 28, and for T = 
and Z — 10. Red-light regions are for AI > 0. The maximum 
value reached is AI ~ 1.4 Blue-dark regions are for AI < 0. The 
minimum value reached is A J ~ —1.4. At the borders of red and 
blue regions the critical currents are equal and AI = 0. 



varying functions of L and M. We here choose the parame- 
ter range 6 < L/c < 28 (g) < M < 0.75 to be consistent 
with both our approximations and with typical experimental 
situations. First, we look at which of the ferromagnetic mech- 
anisms that gives rise to the larger critical current, by consid- 
ering the quantity AI = (/f^^^ - I^'^^)/Iq. This quantity 
is plotted for the chosen range of parameters for T = and 
Z — 10 in Fig. [6] where red-light regions are for AI > 0. 
The maximum value reached is AI ~ 1.4. Blue-dark re- 
gions are for AI < 0, where the minimum value reached 
is AI ~ —1.4. At the borders of red and blue regions, the 
critical currents are equal and AI = 0. The complicated pat- 
tern of AI shows that even a slight change in mid-layer width 
or polarization can reverse the order of critical current mag- 
nitudes in SBAF and STF meaning that AI has a strong local 
character In the low polarization regime both AI < and 
AI > are realized with the same frequency, while in most 
of the intermediate/high polarization regime AI > 0. This 
indicates that a SBAF is likely to support a larger critical cur- 
rent than a STF when considering intermediate/strong ferro- 
magnets. This property is qualitatively independent of width, 
temperature, and interface transparency. 

We now analyze the — tt transitions in the same range by 
plotting phase diagrams for STF and SBAF at T = O.OlTc 
and different values of Z (see Fig. |7]). For the moment, we 
choose to leave aside temperature dependence, since it is not 
the main mechanism driving the — tt transitions.'^ 

From Fig. |7] it is clear that a larger number of M and L 
driven — tt transitions is generally expected for SBAF than 
for STF. In particular, this is always the case for polarization 
driven transitions at fixed width L, regardless of the value of 
Z. For i-driven transitions at fixed polarization, this is the 
case only when M > 0.6. Another point to notice is that 



STF Z=0 SBAF Z=0 




Lk Lk 



FIG. 7: (Color online). M — L phase diagram for STF (left panels) 
and SBAF (right panels) for T = O.OlTc and different Z values. 
White and colored regions correspond to 0— and tt— phases, respec- 
tively. 



for both STF and SBAF the number of possible transitions 
is larger for higher Z values. At this point the origin of the 
complicated patterns seen in Fig. |5] for M and L dependence 
of Ic when Z 7^ is clear since — tt region boundaiies 
become wavy for finite interface barrier strengths such that 
changing only M or L leaving the other fixed, boundaries can 
be crossed at a non-uniform frequency. We underline that the 
phase of the junctions manifests a strong local character and 
can be switched by slightly altering the width and/or the po- 
larization. 

Fig. [T] shows that SBAF undergoes more transitions than 
STF. It is interesting to analyze the number of transitions for 
a general F which can be a mixture of STF and SBAF. Fig. [8] 
shows the number of L driven — tt transitions at fixed po- 
larization for different mixing degree W ranging from STF 
to SBAF. It is shown how the bandwidth asymmetry tends to 
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Number of L driven 0-n transitions 
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W 

FIG. 8: (Color on line). Number of L driven — tt transitions in 
the range 6 < Lk < 28 as a function of the relative contribution of 
SBAF to total polarization. T/Tc = 0.01 and Z = are fixed. 



increase the number of transitions and that this effect is more 
prominent in the high polarization regime where the number 
of transitions is almost doubled. One important point to note 
is that the same number of transitions observed in a SBAF 
(W — 1) can be observed in a mixed F whenever W is large 
enough, e.g. W — 0.5 for M = 0.75. This is an important 
point when relating our results with real materials because in 
experiments, one usually encounters the mixed F case. Even 
if T/Tc = 0.01 and Z = are fixed in Fig. [t] we under- 
line that the general trend of bandwidth asymmetry increasing 
the number of transitions does not depend on their particular 
values. 



IV. CONCLUSIONS 



ization regime. By analyzing Josephson critical current we 
have shown that bandwidth asymmetry can both enhance and 
decrease its value, the former situation being more common 
for strong ferromagnets. Our findings are relevant for many 
interesting magnetic materials which cannot be framed exclu- 
sively within a Stoner scenario. As relevant examples we cite 
the half-metal ferromagnets defined by the property of having 
almost 100% transport spin polarization.'*^ Materials belong- 
ing to this class usually have spin polarization M in the range 
where bandwidth asymmetry is clearly distinguishable from 
Stoner exchange, e.g. Cr02,'^Lao.7Sro.3Mn03,'*'* Fe304,'Sl 
and EuO,^" among others. Their high polarization values can- 
not be ascribed only to Stoner exchange, so that when trying 
to theoretically model them the inclusion of bandwidth asym- 
metry is a fundamental ingredient. 

We would like to emphasize that our findings may turn out 
to be useful in ascertaining if a given ferromagnetic material 
has spin bandwidth asymmetry. Furthermore, since a band- 
width asymmetric F layer can provide more frequent transi- 
tions and is likely to support a larger critical current for a given 
set of parameters than a Stoner ferromagnet, it may be more 
suitable for electronics or spintronics applications and devices 
based on the Josephson effect and relying on — tt transitions. 
Moreover, considering experimental and sample production 
limitations, the choice of magnetization mechanism may work 
as an extra degree of freedom which by itself may determine 
the ground state phase difference across the junction or the 
magnitude of Josephson critical current. 
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In conclusion, we have analyzed the Josephson effect in 
short, ballistic single channel S/F/S junctions taking into ac- 
count the possibility for ferromagnetism to be driven by a 
mass renormalization of carriers with opposite spin, i.e. a spin 
bandwidth asymmetry. We have compared a junction with 
this unconventional kinetically driven ferromagnetism in the 
F layer with one with the usual Stoner mechanism consider- 
ing also their interplay. Analyzing Andreev levels, free energy 
and currents, we have shown that the Josephson effect in the 
two junctions shows different features especially for interme- 
diate/high polarization values. In particular, we have shown 
that the junction with total or partial spin bandwidth asymme- 
try in the F layer undergoes a larger number of 0— tt transitions 
driven by variations injunction width and polarization. By ex- 
amining free energy and phase diagrams, we have pointed out 
how junctions with different magnetic mechanisms in the F 
layer can be in different phases even if all junction parame- 
ters have the same values. We have remarked that this is a 
common (rare) situation in the intermediate/high (low) polar- 



Appendix A: Explicit form of Andreev levels 

The explicit form of A, B, C, D in Eq. (|6]) is 

A = 2{a~b) , 

B = 4:C, 

C==2(a + 6)-d, 
D = 2d , 

where 



(A^ + A|) - 2 



(Al) 
(A2) 
(A3) 
(A4) 



tan (kLAcr) tan (kLAg-) ' 



(A5) 
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4Z^ 



2Z{1 



-2Z xi 



tan {kLA^) tan (fcLA^) ^ ^"^'^ 



A. 



tan(fcLA(j) tan(A:iA5^ 



A„ 



(1 + Z^)( 



' tan {kLAa 



+ A 



A^ 



tan (kLAg-) J ' 



A. 



(A6) 



tan(fcLAcr) tan(fciA5^ 



+ ^^^^ + ^^tan(fcLA.) " ^'^^ + ^^an(fcLA^)' 



and 



where 



d = 8 



X„ 



X-a 



sin (kLAa-) sin (kLAg-) ' 



I m f U 



(A7) 
(A8) 

(A9) 
(AlO) 



Appendix B: Determination of magnetic parameters for mixed F 

Since the polarization is not a separable function with re- 
spect to the exchange interaction and mass mismatch, one 
cannot immediately obtain the values assumed by these mi- 
croscopic parameters when F is a mixture of STF and SBAF. 
In this appendix we illustrate how we perform this calculation 
for a given M and W, where the mixing parameter W quan- 
tifies the relative weight of the two contributions to M. We 
notice that STF coiTesponds to T4^ = and SBAF to = 1. 

Let us consider the two dimensional {m^/mi,U/Ep) 
space parameter over which AI is defined (see e. g. Fig 



1 of RefPS). For a given value of the polarization M, and 
from Eq. ( [TT| , STF is represented by the point 



2M 



,0 



and, from Eq. ( [TO| i, SBAF is represented by the point 
1, (^If^x^ ^ • The distance between these two points along 
a constant magnetization path 1{M) can be written as 



1{M) = 




where B is the incomplete Euler function 



B{z;a,b) = [ t° 
Jo 



Then the coordinates in parameter space associated with a 
mixture defined by the actual W can be evaluated by imposing 
that the path between the STF and the sought points is a frac- 
tion of the total length equals to W 1{AI). Consequently, the 
mass mismatch value for a degree W mixed F at polarization 
AI can be estimated by numerically solving the equation 



„z7r/4 




mi 



AP 



1 



AP - 1 



-1/4,3/2 



; -1/4, 3/2 



Wl{M), 



where 



M2 - 1 



Once the sought value of mass mismatch is known, the cor- 
responding value of exchange interaction can be obtained by 
inversion of Eq. (|9|. 
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